Renormalization approach to constituent quark models of Quarkonium* 
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Constituent quark models, while successful, require a great deal of fine tuning of the short distance 
interactions by introducing phenomenological gluonic form factors which are ultimately designed to 
accurately reproduce the spectrum. We apply and develop renormalization ideas to reduce the 
short distance sensitivity and show that, as naively expected, but not explicitly implemented in 
the models, the physics of binding is entirely linked to the string tension whereas leptonic decays 
depend more on the gluon exchange potential. We also show how the spectrum of S- and D- l - " 
states is successfully intertwined through the singular tensor interaction. 
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I. INTRODUCTION 

After the discovery of the first heavy-quark bound 
states, the ip and the T systems, it was soon realized that 
a nonrelativistic picture seemed to hold for them. The 
charmonium system has become the prototypical 'hydro- 
gen atom' of meson spectroscopy and therefore it was 
the first and simplest case where a bound state of quarks 
could be studied. The old Cornell potential [l[ 



V(r) = 



4a s 
" 3r 



(1) 



provides a rough estimate of the 1 charmonium spec- 
trum and in particular of the J/ip resonance (see e.g. 
Ref . [U, Q for a review and references therein) . In poten- 
tial models, One-Gluon Exchange (OGE) short distance 
singularities appear and phenomenological gluonic form 
factors must be introduced (see e.g. [J] for an early pro- 
posal). This triggers an unpleasant short distance sensi- 
tivity, and as we will show it mainly hides the fact that 
the ground state is actually being used as an input rather 
than a prediction. In this paper we want to make this 
statement more quantitative and to analyze to what ex- 
tent can one disentangle the physics of the ground state 
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to that of the excited states. The method which we will 
be using is based on renormalization ideas. To provide 
a proper and broad perspective and to avoid confusion 
with other approaches closely related, it is of interest to 
review some relevant developments in what follows. 

For the lowest-lying bound states the size of the system 
is smaller or comparable than 1/Aqcd [H, H|, whereas 
excited states start feeling the long range string tension. 
In the case of charmonium this is so until one reaches 
the DD production threshold. Of course, one expects 
that at short distances corrections are computable within 
perturbation theory due to asymptotic freedom, whereas 
at long distances one must resort to lattice calculations. 

For heavy particles the static energy, E(R) — 2m, is 
an observable and thus a gauge and scale independent 
quantity. The rigorous non-perturbative definition of the 
QQ potential is the expectation value of the Wilson loop 
and has been continuously used on lattice calculations (a 
review on potentials from lattice QCD can be found e.g. 
in Ref. Q ) . It should be noticed that potentials are com- 
puted for point-like sources, i.e. structureless quarks. 
The determination of the QQ potential has been carried 
out in perturbation theory up to N 3 LO HQ. Looking 
for a scale independent potential generates the Renor- 
malon problem (for a review see e.g. Ref. (T?jj ) which 
causes serious convergence difficulties. Effective Field 
Theory (EFT) methods disent ang le the hard, soft and 
ultrasoft scales perturbatively [Tl[ and explicitly solve 
the problem. 

The pNRQCD approach 0, Q3] is an EFT which has 
better convergence properties than NRQCD since it en- 
codes the ultrasoft scales non-perturbatively in terms of 
potentials which enter as Wilson coefficients and should 
be obtained directly from the lattice or string models en- 
abling a reliable description of excited states. The bot- 
tomline seems to be that perturbative potentials should 
be treated perturbatively whereas non perturbative po- 
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tentials should be treated non-perturbatively. 

Following an observation 14] that the renormalon 
problem is triggered by an illegal extension of Fourier 
integrals to all momenta, a parametrization of the lattice 
data in the short distance region, r < 0.1 5fm, has been 
proposed [lj| cutting off low momenta below /z/ ~ lGeV. 
If fif —J- a flagrant lack of convergence of perturbation 
theory has been pictorially represented in Ref. [15j; de- 
pending on the scheme prop osed one may even violate 
fundamental inequalities |l6l[l7| . It should be noted that 
this representation smoothly matches the current lattice 
determination of the static energy allowing an analysis 
of the quarkonium spectrum, but suggests that the non- 
perturbative contribution corresponds to q < fif, larger 
than the ultrasoft scales expected from pNRQCD. Un- 
like the fif — > case the convergence is very good and 
basically a simple OGE potential is left. 

Of course all these considerations could less trans- 
parently be handled from a fundamental point of view 
by solving the theory ab initio. However, direct lat- 
tice QCD calculations of quarkonium spectrum require 
small enough lattice spacings, so that Ma -C 1 and so 
far promising calculations operate at a ~ 0.1 — 0.2fm 
producing some over-estimation of excited states [la ]. 

Turning to a less fundamental approach such as the 
non-relativistic quark model and in the spirit of the Cor- 
nell potential and extensions of it let us mention that the 
standard regularization of the short distance contribu- 
tions to the OGE potential (corresponding to delta con- 
tact interactions and the tensor and spin-orbit pieces), 
actually corresponds to tune the parameters of the reg- 
ulator to fit the J ftp mass. In addition, one of the usual 
parameters which are adjusted is the strong coupling 
constant, a s . However, the world average strong cou- 
pling constant is a s {M z o) = 0.1184 ± 0.0007 The 
measured values of a s (Q 2 ), covering energy scales from 
Q = M T = 1.78GeV to 209GeV, exactly follow the en- 
ergy dependence predicted by QCD and therefore signif- 
icantly test the concept of Asymptotic Freedom. There- 
fore, the value of a s at the charm quark mass is fixed and 
is not larger than 0.2 — 0.3. This value contradicts some 
quark models. 

A nonperturbative determination of the QCD potential 
at 0{l/m) was undertaken for the first time in (20l. l2lj|. 
they found a potential which depends on the interquark 
distance as —c/r 2 and it is comparable with the coulom- 
bic term of the static potential when applied to charmo- 
nium and amounts to one-fourth of the coulombic term 
for bottomonium. 

Lattice calculations of 1/ro corrections based on the 
Bethe-Salpeter amplitude [22| (see also (23|) reveal that 
the dependence with the interquark distance of the 1 /to 
contributions to the potential is the same than in the 
case above but other term which grows logarithmically 
with r was found. When the data are fitted taking into 
account LO and NLO contributions to the potential, the 
coefficients which accompany Coulomb and string forces 
are comparable to those calculated through Wilson loops. 



The 0(l/m 2 ) spin-dependent and momentum- 
dependent corrections are investigated in SU(3) lattice 
gauge theory [24j| . These corrections are relevant ingre- 
dients of an effective field theory for heavy quarkonium 
called potential nonrelativistic QCD. 

In the present paper we want to address the question 
on how much can we deduce for the charmonium spec- 
trum from the knowledge of the potential at long dis- 
tances in a way that our ignorance at short distances 
needs not play a crucial role. This allows to disentan- 
gle the physics of the ground state to that of the excited 
states which are sensitive to the perturbative (Coulomb 
like) and non-perturbative (String like) corrections. Ac- 
tually, our original motivation for the present study was 
to analyze the role of regulators within constituent quark 
models like that of Ref. [25[. However, our ideas can be 
generally extended to any quark model and shed some 
light on their predictive power. 

The plan of the paper is as follows. In Section |H] we 
provide a simplified and pedagogical discussion of our 
main points in the simplest case of the Cornell potential. 
The role of singular tensor and spin-orbit interactions is 
analyzed in Section IIIII where an interesting correlation 
between S- and D-waves is found. In Section HVl we re- 
view the bases of the constituent quark model proposed 
in Ref. [H}. In Section [V] we focus on the new renormal- 
ized charmonium model where the short distance regu- 
lators have been removed as they turn out to be phys- 
ically irrelevant for most observables. Interesting issues 
regarding the scope and applicability of these renormal- 
ization ideas are scrutinized in Section IVII Further as- 
pects are dealt with by an enlightening analysis of the 
Bosonic String Model (BSM), Section IVTA"! Finally, we 
give some conclusions in Section IVIII 

II. RENORMALIZATION OF THE CORNELL 
POTENTIAL 

In this section, we provide a comprehensive discus- 
sion of the renormalization approach as applied to 
heavy quark systems and confining potentials within 
the context of non-relativistic quantum mechanics. The 
Schrodinger equation is given by 

- — !— V 2 * + V{r) * = (M -m q -m q )ty, (2) 

2/Xq g 

where fi qq = m q m q / (ra q + m q ) is the reduced q — q mass 
and the normalization condition J d r\^(r)\ = 1 must 
be imposed on the solution. We use the standard Cornell 
potential 

Acv 

V{r)=—£+or t (3) 

for the purposes of illustration. The more elaborated 
potentials, where relativistic effects and spin-dependent 
corrections are added, will be discussed in later Sec- 
tions. It is remarkable that such a simple potential 



3 



not only captures the relevant physics of the prob- 
lem but also is accurately described by lattice calcu- 
lations (see e.g. J2f|) where Eq. ((3} is favoured with 
4a s /3 = 0.25(1), extremely close to the bosonic string 
model a s = n/16 [2?], HH. It is worth mentioning that 
the lattice calculations of q — q potentials correspond to 
using point sources, i.e., elementary quarks until dis- 
tances comparable to the lattice spacing. A smooth 
transition below a ~ 0.2 fm towards OGE has been ob- 
served. Unless otherwise stated we will take a s — 
a = (420 MeV) 2 = 0.1764 GeV 2 , M J/4> = 3096.916 MeV, 
m c = 1200 MeV and r c = 0.3 - 0.01 fm. 

For the spherically symmetric Cornell potential the to- 
tal relative wave function can be factorized in the usual 
fashion, ^(f) = (u(r)/r)Yi m (f), with u(r) the reduced 
wave function and Y/ m (f ) the conventional spherical har- 
monics. 



A. Review of the Renormalization approach 

Let us consider the standard non-relativistic 
Schrodingcr equation for bound states in S- waves 

- ^<(r) + V(r)u n (r) = E n u n (r), (4) 

where u n {r) vanishes at long distances and the energy is 
defined with respect to the q — q threshold, E n = M n — 
m q - rrig. 

Anticipating our discussion we will assume for definite- 
ness a short distance auxiliary cut-off, r c , below which 
the potential vanishes. This cut-off is just a parame- 
ter which will ultimately be removed while keeping some 
physical condition hxed. In our case we will choose to 
fix the ground state energy to the experimental value. 
In the numerical application we will be concerned with 
the residual cut-off dependence of observables induced by 
such a procedure. 

Using the standard trick of multiplying the Eq. (fj| by 
u m (r) and subtracting the similar equation with (n O 
to), we get for two different energies E n ^ E m the orthog- 
onality relation between their bound state wave functions 

u' n (r c )u m (r c ) - u n (r c )u' m (r c ) 

= 2(j,(E n -E m ) u n {r)u m (r) dr. (5) 

J r c 

Note that usually the regularity condition at the origin, 
u n {r c ) = for r c — > 0, is imposed. Whence orthogonality 
of wave functions with different energies holds. However, 
this is not the only solution to the orthogonality require- 
ment. Instead the common boundary condition, we may 
as well take 

u' m (r c ) = u' n {r c ) 
u m (r c ) u n {r c y 

for any two states, meaning that the logarithmic deriva- 
tive at short distances becomes state independent. In 



particular, choosing the ground state as a reference state 
we get the condition 

"o(^) = <(r c ) ^ 
u Q (r c ) u n (r c )' 



How can this logarithmic derivative be determined? If 
we know the energy of the ground state and the potential 
we may integrate from the long distance region inward 
to deduce u' (r c )/uo(r c ). Once this number is known, we 
may use Eq. (JTJ) to integrate out the excited state and the 
corresponding bound state energy can be fixed by requir- 
ing the wave function to vanish at large distances. Thus, 
such a procedure allows to treat the ground state energy, 
Eq, as an independent variable from the potential V(r) 
and still deduce wave functions and the excited spectrum. 
There is of course the question on how to interpret the 
short distance cut-off, r c . In principle one may look for 
stability at scales below the relevant sizes. Actually, by 
varying the cut-off in this region this is a way of assessing 
theoretical uncertainties . However, there are situations 
where the limit r c — > is rather smooth and induces 
moderate changes in observables. 

A good feature of the present approach is that since 
wave functions are matched at short distances the result- 
ing energies are largely independent on the short distance 
behaviour of the potential 2 . 

The previous discussion has been conducted for 
5-waves and regular potentials, i.e. fulfilling 
lim r ^o '1^(01 < 00 • Higher partial waves cannot be 
renormalized in this fashion as short distances are domi- 
nated by the centrifugal barrier. This also has the benefit 
of diminishing the dependence on the short distance po- 
tential since the short distance behaviour of the wave 
function, u(r) ~ r l+1 , is not controlled by the potential. 

The case of singular potentials has also been discussed 
at length (see e.g. Refs. [29l - [3"lj regarding nuclear or 
atomic systems). The relation to momentum space renor- 
malization and the corresponding Lagrangian countert- 
erms is discussed in Ref. 13 21 . The irrelevance of form 
factors is analyzed in Ref. [3l[ [33| . 



1 In a model where phenomenological form factors are imple- 
mented ad hoc, just to prevent singularities, the assessment of 
theoretical errors could be done by choosing all possible regu- 
larization functions, unless the form factor is known from first 
principles. 

2 The previous approach is well documented in the mathematical 
literature as the theory of self-adjoint extensions of hermitian op- 
erators. This endows the Hilbert space with a common domain of 
functions where the solutions of the Schrodinger equation span a 
complete set. The generation of a new scale independently on the 
potential resembles the well known phenomenon of dimensional 
transmutation. 
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B. Analysis of the linear potential 

In order to provide a clear picture of the procedure to 
be carried out in the present work, we discuss the issue 
within the framework of the simplified model where the 
OGE piece is neglected. For S'-wave charmonium states 
we have the simplified problem 



1 



u cc( r ) + ar u cc{r) = Eu C c(r). 



(8) 



As it is well known, the general solution in the inner 
region is a linear combination of regular and irregular 
Airy functions 



u C c(r) = ciAi(z) + c 2 Bi(z), 

where the dimensionless variable 

m c (—M + 2m c + ra) 
(m c a) 2 / 3 



(9) 



(10) 



has been introduced. At large values of the argument one 
has 



Ai(z) = 



Bi(z) 



_2 3/2 

e 3 



2^ 1/4 

2 3/2 



+ 0(z- 3 / 2 ) 



(11) 



so that for bound states and assuming no opening of me- 
son decay channels we may discard the Bi(z) function. 
The standard approach consists of requiring the regular- 
ity condition at the origin. 

Indeed, the regular solution at the origin requires find- 
ing the zeros of the Airy function, Ai(z n ) = 0, yielding 
the quantization formula 



m c (-M n + 2m c ) 
(m c a) 2 / 3 ' 



(12) 



where z = -2.33811, z x = -4.08795, z 2 = 
-5.52056, z 3 = -6.78671, z 4 = -7.94413, ... are the cor- 
responding lowest zeros. Note that in this particular 
case the string tension and the quark mass determine 
the binding energy completely. 

An alternative procedure might be as follows. Since 
we want to fit the J/tjj mass and we are confident that 
the long distance dynamics is given by the linear poten- 
tial Q, we are quite certain that at long distances the 
wave function is given by 



; (r) = ciAi(z), 



(13) 



where now Mj/^ is taken as an input. Of course, for a 
generic value of the string tension the wave function will 
fail to vanish at the origin (see Fig.[T]). This may be seen 
as a drawback but the gain is in the prediction of the 
excited spectrum. This is done by matching logarithmic 
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FIG. 1: Wave functions for the different states. The normal- 
ization is such that at the matching point, r c — > Ofm, the 
functions and derivatives of the different states coincide. 



derivatives at the origin or at a short distance cut-off, 
Eq. 0, yielding 



Ai'(zo) _ Ai'(* n ) 



M(zo) 
where z n is given by 



Ai(z„) ' 



m c (—M„ + 2m c + ar c ) 



{ m 



c a)V3 



(14) 



(15) 



This condition guarantees the orthogonality of states 
and provides M n from Mq for any value of r c . Note 
that by definition the short distance behaviour of the 
wave functions is very similar, so that even though we 
may not know accurately the potential at short distances 
there is an increasingly large cancellation. The cut-off 
dependence of the states is depicted in Fig.[5] 3 . As we see, 
changing the scale provides a mild dependence featuring 
the advertised short distance insensitivity. 



C. Inclusion of One-Gluon Exchange 

In heavy quark systems the OGE potential is consid- 
ered to be a very important ingredient to provide some 
additional short range corrections. Here we will show 
that as far as binding is concerned this OGE contri- 
bution becomes almost irrelevant. Contrary to what 
one might naively think this is actually good news. In- 
deed, as we mentioned in the introduction the short dis- 
tance part of the potential is not well known even at 



3 The cut-off dependence can be understood on purely analytical 
grounds as carried out in Appendix [A] for a wide class of poten- 
tials. 
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FIG. 2: Dependence of the excited charmonium states on the 
short distance cut-off. 

the OGE level. If we write the OGE potential in the 
form V(r) = — 4a s (/i)/3r, with /i the MS renormaliza- 
tion scale, it is not obvious what scale should one use a 
priori. Lattice calculations display such a behaviour at 
short distances (see the recent update where running is 
actually observed [34[ ) . There have been attempts to im- 
prove on this by imposing renormalization group invari- 
ance of the potential. Unfortunately such a procedure 
breaks down due to the appearance of renormalons (see 
e.g. Ref. [To| and [lj|) which spoil a convergence pattern. 
In addition, there is also the problem on how far should 
this OGE force be extended, since it only applies to very 
short distances. However, one expects that as compared 
to the linear potential the long range effect is not crucial. 

The region where OGE becomes comparable with the 
linear potential is r ~ b = y ot s /a ~ 0.2 — 0.3fm where 
the total potential is about V(r) ~ — 180MeV. At much 
larger distances we may neglect the OGE component so 
that we may take the asymptotic behaviour of Eq. (|13[) . 
The resulting wave functions without and with OGE arc 
depicted in Fig. El Similarly to the case without OGE the 
procedure is as follows. We start from the given ground 
state energy Mq = Mj/^ integrating inward down to the 
short distance cut-off radius, r c , the full Cornell poten- 
tial, Eq. ©. To facilitate the comparison we take the 
same asymptotic wave function for both cases, without 
or with OGE. As can be seen, the OGE attraction pro- 
vides a larger probability in the inner region below the 
scale where OGE starts being negligible. Our numeri- 
cal calculations indicate that this long distance scale is 
about 1.5fm. From the ground state we can construct the 
remaining excited states by matching logarithmic deriva- 
tives at the short distance cut-off, as indicated above. 

The short distance cut-off dependence of the excited 
charmonium states for several values of the strong cou- 
pling constant is similar in all cases to the a s = situ- 
ation shown in Fig. [5J While at large r c the calculation 
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FIG. 3: Charmonium wave functions for ip(lS) and ij){2S) 
states without (solid line) and with (dashed line) One Gluon 
Exchange potential. 
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FIG. 4: Mass shift (in MeV) as a function of the strong cou- 
pling constant a a when the short distance cut-off is removed 
r c -> for the ip(2S), ip(3S) and tpi^S) relative to the ground 
tp(lS) state. 



is dominated by the linear potential, wave functions at 
short distances are very much alike due to the common 
boundary condition. Thus, the difference in energy due 
to OGE comes from the energy dependence of the wave 
function below the non-perturbative scale. In Table U 
for the case r c — > 0, we show numerical values of the ex- 
cited states for a s = 0.0, 0.1, 0.2 and 0.3. The mass shift 
AM n = M n - M is shown in Fig. H for r c = 0.06fm and 
as we can see the OGE effects are rather moderate. 
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3615 


3630 


3638 


3640 
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4039 


4060 


4070 


4073 


3 


4414 


4439 


4449 


4452 


4 


4756 


4783 


4795 


4798 



TABLE I: Masses in MeV of S-wave charmonium states for 
several values of the a s coupling constant and keeping always 
the ground state mass M = 3096.916 MeV. 



III. COUPLED CHANNEL SOLUTIONS: 
MIXING S AND D STATES 



In this section we include the spin dependent potential 
contributions coming from the OGE interaction to the 
simple Cornell potential discussed in Section HT1 



The radial Schrodinger equation for the 3 Si — 3 Di cou- 
pled channel reads 

- u"{r) + U s (r)u(r) + U SD (r)w(r) = n 2 u(r), 

6 



w"(r) + U S D{r)u{r) + 



U D {r) 



w(r) — k w(r), 
(20) 



where k 2 = m c (M — 2m c ). Us, Ud and Usd are the 
different contributions of the reduced potential, U(r) — 
2/j,V(r), where the 3 Si — 3 Di coupled channel potential 
is given by 



V s (r) = 
V D (r) = 

V SD (r) 



4as 
2y/2 a s 



ar, 



20 a s 1 



(21) 



A. Potential 

The qq interaction can be written as 
V{r) = ar + V GE{r), 



(16) 



where now Voge(^ contains central, the coulomb inter- 
action, and non-central, tensor and spin-orbit, contribu- 
tions 



VP 



r S G E(f) 



4a s 1 
3 r' 
1 «» 1 

— ■ S . 

m r° 



(17) 
(18) 
(19) 



where a s is the strong coupling constant, m is the equal 
quark and antiquark mass, = 3 (t?i •f)(tT2 ■f) — &x-&2 is 
the quark tensor operator, S the total spin operator and 
L is the relative orbital angular momentum operator. 



B. The coupled S-D equations 

Our previous discussion of renormalization was under- 
taken without taking into account the role played by the 
tensor and spin-orbit forces. In the case of charmonium, 
these states are a combination of S and D-wave compo- 
nents due to the tensor force contribution. As we will 
see below, this tensor force is small enough to have al- 
most pure S and D orbital state components. The inter- 
esting aspect of our discussion below is that using just 
one renormalization condition we can predict all S and 
D-wave mesons, i.e. we reduce the number of renormal- 
ization conditions. 



in which a s is the strong coupling constant, a is the string 
tension and m c is the quark mass. Obviously, in order to 
describe a bound state we seek for normalizable solutions 



[u{rf +w{r) 2 ] = 1, 



(22) 



which impose conditions on the wave functions both at 
infinity as well as at the origin. 

The set of equations ([20| must be accompanied by 
asymptotic conditions at infinity. As in previous sec- 
tions, we have a linear potential at large distances for 
both channels and therefore the wave function is a lin- 
ear combination of regular and irregular Airy functions. 
Once we have discard the irregular function at long dis- 
tances, Bi(z), the wave functions at infinity have the fol- 
lowing behaviour 



u(r) -> A s Ai(z), 
w(r) -> A D M{z), 



(23) 



with z defined in Eq. (jlOl) and As is the normaliza- 
tion factor and the asymptotic D/S ratio parameter is 
defined by rj — An I As- Ideally, one would integrate 
the Schrodinger equation taking its solutions at infinity, 
Eq. (|23[) . which depend on the bound energy and rj. The 
singular structure of the problem at short distances re- 
quires a specific analysis of the cou pled equations as it 
has been done extensively elsewhere [33j and we adapt in 
Appendix [C] for our particular case. The result amounts 
to integrate from infinity for the physical value of Mj/^ 
and rj. Generally, the solutions diverge strongly at the 
origin, so that the normalization of the state is precluded. 
However, there is a particular value of rj which guaran- 
tees that the wave function becomes normalizable. Then, 
if one imposes the regularity condition at the origin one 
will determine rj and therefore the wave function of the 
bound state. In practice, however, the converging solu- 
tion is rather elusive since integrated-in solutions quickly 
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run into diverging solution due to the round-off errors 
and dominate over the converging solution. 

According to Ref. [33| one may proceed as follows. One 
can impose different auxiliary short distance boundary 
conditions corresponding to a choice of regular solutions 
at the origin 



u(r c 
u'(r L 
w(r c 
w'(r c 

u(r c ) — v / 2w(r c 
u'{r c ) - V2w'(r c 



= (BC1), 

= (BC2), 

= (BC3), 

= (BC4), 

= (BC5), 

= (BC6). 



(24) 



All these boundary conditions must predict the same 
value of 77 at some value of the cut-off radius, r c , the 
precise convergence value corresponds to the particular 
choice. As in other Refs. [35| we find convergence for 
the boundary conditions BC5 and BC6 for larger cut-off 
radii, which improves the numerical results. 

To calculate the D/S asymptotic ratio, 77, it is con- 
venient to use the superposition principle of boundary 
conditions [35] to write 



u(r) = u s (r) + rju D (r), 
w(r) = w s (r) + T]w D (r), 



(25) 



where (us, u>s) and (u_d, wd) correspond to the boundary 
conditions at infinity, Eq. ([23]) with As = 1 and An = 
and with As = and Ajj = 1, respectively. Through 
that decomposition the boundary conditions, BC1-BC6, 
can be rewritten as algebraic expressions for 77. For in- 
stance, if we use the BC6 boundary condition we get 



n 



u' s (r c ) - \/2w' s (r c ) 
i' D {r c ) - y/2w' D (r c ) 



(26) 



Once t\ has been calculated, the wave function of the 
bound state is completely determined by the normaliza- 
tion factor As 



u(r) = A s (u s +r]u D ), 
w(r) = A s (w s + V w d), 



(27) 



in which As is obtained normalizing the wave function 
to one 



A s I [{u s + riu D f + (w s + riw D ) 2 ]dr = 1. (28) 



The above procedure can be undertaken for the ground 
state of the system if its energy is known. Now, if we 
want to calculate the excited states of the system we 
must impose the orthogonality condition between wave 
functions of states with different energy together with 
the regularity condition at the origin. 

For a regular potential the orthogonality between wave 
functions of states with different energy is a property of 



the solutions of the Schrodinger equation and it always 
holds. However, if the potential contains some attractive 
singular contribution we have to impose explicitly the 
orthogonality between wave functions. 

Thus, given the ground state and one excited state, the 
orthogonality condition can be written as 



dr [u Q (r)u m (r) + w Q (r)w m (r)} = 0, (29) 



where it is useful to rewrite the above expression through 
a Lagrange identity 



[u' u m - u u' m + w' w m - w w' m ] 



= 0. 



(30) 



Note that any individual term in the integrand is actu- 
ally divergent, because of the dominance of the singular 
solutions at the origin. At very short distances, the or- 
thogonality between wave functions and the regularity 
condition of them have been imposed at a certain cut-off 
radius, r c . Of course, we always check that the numerical 
calculation is stable against suitable changes of the short 
distance cut-off so that the range r c ~ 0.01 — 0.3fm is 
sufficient. In that situation, the orthogonality condition, 
Eq. (|3"0")l , can be written as 



u' (r c )u m (r c ) + w' (r c )w rn (r c ) 
= u {r c )u' m (r c ) + w (r c )w' m (r c ), 



(31) 



and combining this expression with the corresponding 
one of the boundary conditions, Eq. (|24p . we obtain in 
the case of the boundary condition BC6 



w' {r c ) 



V2u m (r c ) + w m (r c ) y/2 u (r c ) + w (r c ) 



(32) 



and similarly for all other auxiliary boundary conditions. 
Obviously in this case the D/S mixing of the excited 
state is determined from the requirement of regularity at 
the origin, 



u D, m ( r c) ~ V2w' D ^ m (r c ) 



(33) 



We will also compute other observables. The results 
are presented in Table [TTJ In the numerical integration a 
certain maximum radious has to be selected to imposed 
the asymptotic boundary conditions. Due to the long 
range coulomb force the values of the asymptotic param- 
eters are slowly convergent and we quote the values for 
f max = 20 fm. The corresponding wave functions can be 
looked up in Fig. O 

As we see the description of the spectrum is not partic- 
ularly accurate, but one should take into account that we 
are only providing the J ftp mass as an input, as well as 
a s , m c and a, for which we take generally accepted val- 
ues. The present analysis suggests that exploiting these 
ideas one may find a reduction of parameters in quark 
models. As we will see in later Sections, rather than 
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being a mathematical curiosity, these correlations which 
are unveiled by the renormalization approach actually 
are embodied in quark models in a much less transpar- 
ent way. 



IV. CHARMONIUM WITH FORM FACTORS 



Regarding the confinement interaction 4 , A is a global 
constant adjusted to fit the origin of energies, a c and fi c 
are model parameters. At short distances this potential 
presents a linear behavior with an effective confinement 
strength a = —a c /i c {X\ -X^) while it becomes constant at 
large distances. This type of potential shows a threshold 
defined by 



In the previous sections we have presented how one 
may incorporate a desirable short distance insensitivity 
into several simplified models for quarkonium. In this 
section we give a brief description of the specific char- 
monium model which fits a wider and successful phe- 
nomenology (see Ref. [25[ for further details) paying spe- 
cial attention to the form factors introduced to regulate 
the unpleasent short distance singularities. Our intention 
is to re-analyse the model by trading these form factors 
into a less model dependent set of renormalization condi- 
tions. While the model furnishes unequal quark species 
we will restrict for illustration purposes to the case of 
charmonium. 

The corresponding potentials for the qq system stem 
from the non-relativistic reduction of OGE and the con- 
finement component. We separate the central, tensor and 
spin-orbit pieces, as follows: 



One-Gluon Exchange 



V8ge{?) 



1 



V, 



1 a 



6 m 



I e -r/r {n) 

(<?1 • <7 2 )■ 



rrg(/i) 



OGE 



« = -16^'^ 



1 e ~ r / r s M 



ySO 
v OGE 



(r) 



r 

3 a 
8m 
1 



3r 2 (p) rr g (n) 



Sl2, 



2 (Af - A|) 

-r/r g (n) 



L-S. 



• Confinement 



V, 



CON 



vS8n{^ = -{\i-\ c 2 ) 



(f) = -«c(l-e-^) + A (Af-A§) 



a c ^ c e 



Am 2 i 



-(6-8a s )L-S, 



(34) 



(35) 



where ro(fi) = ro^^jf- and r s (/i) = fg 1 ^- are short dis- 
tance regulators and (i stand for the reduced quark mass. 
The contact term of the central potential of one gluon ex- 
change is regularized in a suitable way by replacing the 
Dirac delta function by a Yukawa form 



S(f) 



— > 



1 



,—r/ro 



Airr 2 



(36) 



In table IIIII we show the model parameters constrained 
by the light quark phenomenology and also appear in the 
potentials operating in the heavy quark sector. 



V thr = {~a c + A}(\1 ■ X c 2 ). 



(37) 



By construction, no qq bound states can be found for 
energies higher than this threshold. Actually, the po- 
tential suffers a transition from a color string configura- 
tion between two static color sources into a pair of static 
mesons due to the breaking of the color flux-tube and 
the most favoured subsequent decay into hadrons. In a 
more general set up, a dynamical coupling to the formed 
meson pair should be included. 

The wide energy range covered by a consistent descrip- 
tion of light, strange and heavy mesons requires an effec- 
tive scale-dependent strong coupling constant that can- 
not be obtained from the usual one-loop expression of the 
running coupling constant which diverges at Q 2 = kq CD . 
Following previous work (25|, we use a frozen coupling 
constant such as 



hi 



(38) 



where /i is the reduced mass of the qq pair and ao, Ho and 
Ao are parameters of the model determined by a global 
fit to the hyperfine splitting from the light to the heavy 
quark sector. In our case this yields a value of a s = 0.288. 

We note that the regulators in the tensor and spin- 
orbit terms of the one gluon exchange potential as well 
as the smeared (5-function of the central component are 
introduced just because tensor and spin-orbit terms be- 
come singular at short distances. As a consequence the 
collapse of the system which would occur for the pure un- 
regularized ^-function is prevented. The previous model 
describes successfully the meson spectroscopy [25[ and 
provides a good agreement with the experimental data 
of different physical observables as the strong and radia- 
tive decays of charmonium performed in Ref. [371 ] . 

As we have discussed above it is interesting to inquire 
the real need of form factors which, while providing an ac- 
ceptable phenomenology, are introduced ad hoc to avoid 
the short distance divergences. Indeed, if the short dis- 
tance regulators are removed, i.e. rQ,r g — » 0, the ex- 
pressions of our potentials become Eq. (TIT)]) for the OGE 



4 Our expressions imply the validity of Casimir scaling both for 
the perturbative (tree level) as well as the non-perturbative (con- 
finement) contributions to the potential. A recent perturbative 
calculation to three loops displays violations to this Casimir scal- 
ing |36| . although nothing is implied for the non-perturbative 
contribution. 
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State 


M (MpV\ 

1\± \ 1VJ.C V J 


'1 /is 


(fm- 1 / 2 ) 




P n (%) 




3096.916+ 


-0.0037 


4.01 


0.637 


1.52 


i>(lD) 


3577.7 


-19,71 


0.084 


0.93 


99.91 


4>(2S) 


3634.2 


+0.015 


3.05 


1.05 


1.24 


i/>(2D) 


3995.9 


-17.73 


0.104 


1.29 


99.84 


i>(3S) 


4065.1 


+0.028 


2.66 


1.40 


1.15 


i>(3D) 


4368.7 


-15.70 


0.12 


1.60 


99.75 


i>(4S) 


4443.1 


+0.039 


2.44 


1.70 


1.12 


V>(4£>) 


4710.8 


-14.09 


0.13 


1.88 


99.65 


4>(5S) 


4787.4 


+0.048 


2.29 


1.98 


1.14 



TABLE II: Different observables computed through the coupled toy model and applied to charmonium spectrum. 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 3.5 4 

r(fm) rfftn) r(fm) r(fm) 



FIG. 5: S and D-wave functions of the different states of charmonium calculated through the coupled toy model. 



Quark mass 


m c (MeV) 


1763 


Confinement 


a c (MeV) 


507.4 




li c (fm -1 ) 


0.576 




A (MeV) 


184.432 




a s 


0.81 


OGE 


ao 


2.118 




A,, (fm" 1 ) 


0.113 




fio (MeV) 


36.976 




fo (fm) 


0.181 




f g (fm) 


0.259 



TABLE III: Quark model parameters. 



piece and the confinement part remains the same as in 
Eq. (1551) as this latter contribution does not depend on 



regulators. We will use these expressions for the poten- 
tials in the renormalization scheme. Note that we have 
also discarded the Dirac delta functions. While this may 
seem weird, these are distributions which are not seen in 
the compact support test functions implicitly implied by 
the boundary condition regularization. This also applies 
to any derivatives of the Dirac delta function. This re- 
sult was suggested [3(| and explicitly checked by using a 
momentum space regularization with so-called countert- 
erms 1321] . 



Moreover, the form factors are naturally adjusted in 
order to reproduce some well established physical observ- 
ables, such as the J/ip mass and its leptonic width. What 
will be shown below is that one can actually treat these 
observables as suitable renormalization conditions, with- 
out any specific need of form factor regulators. 
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r c (fm) 




FIG. 6: Dependence of the mass (in MeV) 2 3 Si (upper panel) 
and 2 3 Z?i (lower panel) excited charmonium states on the 
short distance cut-off, r c (in fm). 



V. RENORMALIZED CHARMONIUM MODEL 
A. Masses (uncoupled case) 

In this section we carry out the renormalization, pre- 
sented previously along the lines described in detail in 
Section |nl for the unregularized model potential. While 
the tensor force induces a mixing between S and D- waves 
we will in a first step neglect such a mixing. This implies 
that both ground states are completely unrelated and 
renormalization is pusued independently. As expected, 
the charmonium masses in this scheme depend on the 
short distance cut-off, r c . Fig. [6] shows this dependence 
for the first radial excitations of S and D-wave charmo- 
nium states. One can see that at some value of r c the 
masses do not depend on the short distance cut-off. 

In Tabic IIVI the masses predicted by renormalized 
charmonium model and the standard constituent quark 
model (with form factors) are displayed. We find a per- 



State n 


Mrsc (MeV) M C qm (MeV) 


M exp (MeV) 


Ref. 


3 Si 1 




3096 


3096.916 ± 0.011 


[38] 


2 


3703 


3703 


3686.093 ± 0.034 


[38] 


3 


4097 


4097 


4U39.6 ± 4.o 


fool 

138] 


i 

4 


4389 


4389 


4oot ± 9 ± y 


139] 


- 




A Q~\ A 

46t4 


a a i a 
46t4 


4634_ 7 _ g 


140] 


3 D 1 1 


3796+ 


3796 


3772.92 ± 0.35 




2 


4153 


4153 


4153 ±3 


m 


3 


4426 


4426 


4421 ± 4 


[38] 


4 


4641 


4641 


4664 ±11 ±5 


[39] 



TABLE IV: Predicted charmonium masses (in MeV) within 
the renormalization scheme (RSC) and potential model with 
form factors (CQM). 



feet agreement between both schemes. This provides con- 
fidence on the way the original model took into account 
the unknown short distance dynamics, on the one hand, 
and also on the irrelevance of those form factors for ex- 
cited states as long as the ground state mass is fixed. 

As it becomes clear from the expressions of the po- 
tential, the perturbative and non-perturbative contribu- 
tions dominate at short and large distances respectively. 
Therefore, we want to study now the dependence of the 
mass with respect to two important model parameters, 
the strong coupling constant, a s , and our effective string 
tension of the confinement potential, c, featuring those 
short and long distance effects. 

The dependence on a s while fixing the threshold for 
confinement, Vthr in Eq. (f3"T|). can be seen in Fig. [7] for S 
and D-wave states. For S'-wave states the mass changes 
about 15 MeV meaning a 0.4% in the total mass and 
we see a flattening behaviour along the considered range 
of a s . We find a similar trend for the D-wave states, 
although in this case the change in mass is larger than 
in the case of S'-wave states, the picture shows again a 
rather clear plateau. 

We now turn to the mass dependence on the effective 
string tension of our confinement potential. Again and 
for clarity of presentation we fix the threshold for confine- 
ment, Vthr in Eq. (|37|) . Fig. [8] shows such a dependence 
for the S and D-wave states. The range for the effective 
string tension is in percentage level equal than the range 
of the strong coupling constant and we can see that the 
masses change on the hundreds of MeV. We can conclude 
that at least the masses of excited states are dominated 
by the confinement potential as long as the ground state 
mass is kept to its physical value. In Appendix [D] we 
further analyze the sensitivity of the renormalized model 
parameters. 
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FIG. 7: Mass (in MeV) dependence on the strong coupling 
constant of excited charmonium states. 



B. Leptonic widths 

We now focus on the leptonic widths. This will 
illustrate the interesting subject of the wave func- 
tion renormalization. A complete calculation of 
^(vector meson) — > e + e~ widths involves radiative and 
rclativistic contributions. In the heavy quark limit one 
can write the leptonic decay for iS-wave states as [4l[ 



T (n 3 S*i 



e+e - ) = ^el\R nS m 7 



M 2 



1 - 



16a< 
~3T 



(39) 

where e c = 2/3 and M n is the mass of the charmonium 
state. Decay widths depend on the value of the wave 
function at short distances. Of course while the decay 
may be triggered by a short distance operator, we may 
predict decay ratios as 



r (n 3 ^ -> e +e~) \R nS (0)\ 2 M* 
r(l3SWe+e-) |i? 15 (0)| 2 M2' 



(40) 



FIG. 8: Dependence of the mass (in MeV) excited charmo- 
nium states on the effective string tension of our confinement 
potential, a — —a c fj, c (Af • AJ) (in GeV 2 ). 



This corresponds in practice to implement a common 
wave function renormalization which factors out in the 
ratio. 

Fig. |ni[a) shows the dependence of 1Z with respect to 
the short distance cut-off for the first radial excitation of 
S'-wave states. At some range of r c , the ratio does not 
depend on its value. Fig. [9f b) shows the dependence of 
the ratio for the S'-wave states along the range of strong 
coupling constant. This range is the same as in the case of 
studying masses. We find a stronger dependence on the 
strong coupling constant as expected since the leptonic 
decay is a short range observable. 



C. Inclusion of coupled channels 

Finally, we take into account the role played by the 
tensor force which, as we know from previous sections, 
is the responsible of the mixing between S and D-waves 
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FIG. 9: Top picture shows the dependence of the ratio 1Z 
with respect to short distance cut-off, r c . Bottom picture 
shows the dependence of the ratio 1Z for the excited S-wave 
charmonium states on the strong coupling constant, a s . 



in the case of charmonium. The interesting aspect of 
applying renormalization procedure to S-D coupled case 
is that using just one renormalization condition, the mass 
of the J/i/j, we can predict all S and D-wave mesons. 

The procedure is the same than in Sec. IIIII with only 
one difference, where our confining interaction is given 
by Eq. (f3"5l) . which flattens out at large distances. So in 
this case the asymptotic boundary conditions become 



u(r) -> A s e Kr , 
w(r) -> A 



D 



3 

KT 



(41) 



with k 2 = m c (Vth. r +2m c — M) and Vthr given by Eq. (pJT|) . 

Results on the mass, the asymptotic D/S ratio parame- 
ter and D-wave probability are presented in Table IVl The 
comparison between renormalization scheme and con- 
stituent quark model with form factors is given in Ta- 
ble I VII One can see that the agreement is completely 



State 


M (MeVI 

J.VJ. I Lv J. C V 1 


'1 


Pn (%) 


J/lp 


3095.78 1 " 


-0.0115 


0.146 


4>(2S) 


3703.19 


+0.0105 


0.221 


V>(LD) 


3795.93 


-12.47 


99.852 


4>(3S) 


4097.46 


+0.0321 


0.446 




4152.20 


-9.49 


99.605 




4389.10 


+0.0522 


0.758 




4425.43 


-7.76 


99.280 




4614.28 


+0.0698 


1.108 


i/>(4L>) 


4640.00 


-6.70 


98.925 



TABLE V: Masses (in MeV), the asymptotic D/S ratio pa- 
rameter and D-wave probabilities of charmonium states in- 
cluding S-D mixture. The only input is Mj/^,. 



satisfactory. Essentially, this proves that the form fac- 
tors only provide the correct mass of J/ifi. Once this is 
fixed the rest of the excited states with either S- or D- 
wave character are predicted. This is the main result of 
the present study. 



VI. ON THE APPLICABILITY OF THE 
RENORMALIZATION PROCEDURE 

In most examples shown in this paper, we have always 
dealt with a situation where it was possible to fix the en- 
ergy of a bound state, so that there is a mild dependence 
of excited states on the short distance cut-off. This scal- 
ing analysis is carried out in Appendix [5] providing an 
understanding of the features found numerically. 

Of course, there is the question when the present renor- 
malization procedure is guaranteed to work, i.e. can we 
always mock up our ignorance at short distances by a 
given renormalization condition? On the other hand it is 
also important to know if the answer is affirmative what 
does the renormalization tell about the true solution? 

It is instructive to analyze first the case where renor- 
malization can be carried out at any level of approxima- 
tion of the "true" potential. We provide some insight 
into the renormalization problem by analy sing in some 
detail the bosonic string model (BSM) [271 128||. Such a 
model provides an alternative approach to the physics of 
confinement from the point of view of large distances (for 
a review see e.g. [42[ )• It should be noted that with the 
constraint a s — 7r/16, the Cornell potential corresponds 
to the large distance expansion of the following potential 



V(r 



7 L -Ti + --- ( 42 ) 
2r 8r 6 



where ur\ = ir/6. With the standard value for the 
string tension ^fa — 0.420 GeV one has r ~ 0.3 fm. A 
feature of this potential is that for r < ro it becomes 
purely imaginary. Moreover, lattice calculations [26[ 
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Meson n 



Renormalized scheme 



Form factors scheme 



Exp. data 



J/i> i 

ip(2S) 2 

V>(3770) 3 

V>(4040) 4 

^(4160) 5 

X(4360) 6 

V>(4415) 7 

X(4630) 8 

X(4660) 9 



V 3Sl (%) V 3Dl (%) M (MeV) 



P 3Sl (%) V 3Dl (%) M (MeV) 



99.85 
99.78 

0.15 
99.55 

0.39 
99.24 

0.72 
98.89 

1.08 



0.15 

0.22 
99.85 

0.45 
99.61 

0.76 
99.28 

1.11 
98.92 



3096 1 " 
3703 
3796 
4098 
4152 
4389 
4426 
4614 
4640 



99.96 
99.96 

0.03 
99.94 

0.06 
99.91 

0.09 
99.88 

0.11 



0.04 

0.04 
99.97 

0.06 
99.94 

0.09 
99.91 

0.12 
99.89 



3096 
3703 
3796 
4097 
4153 
4389 
4426 
4614 
4641 



M (MeV) 

3096.916 ± 0.011 
3686.093 ± 0.034 
3772.92 ± 0.35 
4039.6 ± 4.3 

4153 ±3 
4361 ± 9 ± 9 
4421 ± 4 
4634t^ 
4664 ± 11 ± 5 



TABLE VI: Comparison of different properties of charmonium between the renormalization scheme and the constituent quark 
model with form factors, considering coupled channels in both cases. 



yield 4a s /3 = 0.25(1), extremely close to the BSM choice 
of a s — 7r/16 when the Cornell potential, Eq. is 
taken. It has been found [43[ that the transition from 
perturbative to string behaviour is evident from the data 
and takes place at surprisingly small distances 5 . 

The BSM potential satisfies exact fundamental in- 
equalities for the Wilson loop exploiting reflection posi- 
tivity [H[i3, namely 

V'(r) > V"(r) < 0. (43) 

The latter identity is saturated by the linear potential. 
This is satisfied by any inverse power series with negative 
coefficients, i.e. 

oo 

V(r) = or + V + Y / c n /r n , ( 44 ) 

n=l 

where one has c n < 0. However, note that the converse 
need not be true, namely Eq. (|43]) does not necessarily 
imply that all coefficients in Eq. (|4"4")) are negative. 



A. Renormalization of the bosonic string potential 

In this section we deal with the renormalization of the 
bosonic-string potential given by Eq. (1421) . Of course, this 
potential is questionable below the critical radius where 
it becomes imaginary. On the other hand, if the large 
distance expansion is truncated to a finite order we have 
to deal with increasingly singular power divergences. We 
will illustrate how the renormalization procedure works 



5 This of course raises the question on whether the perturbative 
OGE contribution,— 4a s (fi)/3r, should be added on top of the 
BSM-l/r result or it is actually dual to it, in which case a scale 
of p, — 2GeV is implied, if one demands a s (p) = tt/16 at LO in 
perturbation theory. 



when the full potential is compared to the series of sin- 
gular potentials at the origin. 

This example also shows how to treat singular poten- 
tials. Indeed, at NNLO one encounters a 1/r 3 short dis- 
tance divergence. The problem with such a singularity 
is that the energy becomes unbound from below. To 
overcome this situation one can fix, as it was discussed 
in Sec. [TTl the ground state energy and then, using the 
orthogonality condition, the excited states can be pre- 
dicted. Of course, this procedure can also be carried out 
without a singular potential, but just on the basis that 
the short distance behaviour of the potential may not be 
fully reliable. 

Our numerical results are displayed in Fig. [TU] There 
we show the BSM potential compared to the linear 
VLo(r) = or (LO), V Nh0 (r) = or + c x /r (NLO) as well 
as VnnloM = or + cx/r + c^/r 3 (NNLO) approxima- 
tions. Actually, for the typical value of the string tension 
the difference between the potentials is at most about 
50MeV. Applying the renormalization procedure we see 
that when the short distance cut-off coincides with the 
critical distance the difference among the different masses 
become much smaller. This is just a manifestation of the 
short distance insensitivity embodied by the renormaliza- 
tion. While in the exact BSM we are forced to stop at the 
critical radius, the LO, NLO, NNLO, ... truncations to 
such a potential allow to remove the cut-off completely. 
However, as we see from the figure this provides in our 
view an estimate on the systematic error inherent to the 
BSM model; the mass shift generated by the different 
truncations is about 50MeV. Of course, there is no point 
in going beyond such an accuracy before other impor- 
tant effects (like e.g. relativistic corrections) are taken 
into account. 

Actually, there are hints 0, EH that for the Nambu- 
Goto string the critical radius moves about a factor of two 
down (in our case it means ~ 0.2fm ) when quark masses 
become comparable to the string tension m q /y/o ~ 1 (in 
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FIG. 10: Left panel: The bosonic-string model potential (in GeV) and different terms contributing to the long distance 
expansion truncated to a given order as a function of the distance r (in fm). Middle panel: Short distance cut-off r c (in fm) 
dependence of the S-wave charmonium spectrum masses (in GeV) for the different approximations to the potential. Right 
panel: Relative n s Si — > e + e~ decay ratios with respect the leptonic width of the ground state, J/ip —¥ e + e~ . The bands 
indicate the corresponding experimental number. In all cases the vertical line marks the critical radius below which the full 
potential becomes imaginary. Solid line: Full BSM. Thick dotted: LO approximation. Dotted: NLO. Dotted-dashed: NNLO. 



our case m c /y/a = 2.85) 6 

We depict the relative n 3 Si — > e + e~ decay ratios with 
respect the ground state J/tp — > e + e~ decay. As we have 
emphasized above the advantage of such a comparison is 
that purely perturbative and state independent correc- 
tions in the van Royen-Weisskopf formula are factorized 
out. However, note that the very use of this formula as- 
sumes a regular potential. Thus we take naturally the 
wave function at the shortest cut-off radius, r c . As we 
can see there are somewhat large differences for the sim- 
ple linear potential (LO) approximation as compared to 
NLO, NNLO or the full BSM results which agree among 
them, displaying a convergent pattern at any value of r c . 
Differences between NLO and NNLO start being notice- 
able below cut-off radius of about 0.15fm for the first ex- 
cited state, although the value seems slightly displaced to 
larger numbers for the next excited states, as one might 
expect from their stronger short distance sensitivity. Also 
interesting is the fully converged result at r c — > 0, a fact 



6 We also note a recent effective string theory approach [46| | , where 
these finite mass effects are evaluated regardless on the specific 
string dynamics, yielding a potential, up to a constant term, of 
the form 



V(r) = - 



■ log(v / K»") 



(45) 



hence they obtain the nontrivial result that the l/m potential 
must grow logarithmically at large r and it has to connect with 
the 1/r 2 behaviour at short distances of the l/m contribution 
calculated by perturbation theory. The comparison of the lattice 
data is done through the next potential form 



Vfl t (r) 



1 



(46) 



where once a fit to the lattice data is done [20 . |21| | yields c' 
0.0027GeV 2 and d' = 0.0075GeV 2 . 



which follows analytically by a detailed study of the short 
distance wave functions (see Appendix [Cj . 

These features reinforce the value of renormalization 
as a way to study the effect of improving on the long 
distance components of the qq potential. Indeed, even 
though the NNLO potential strongly diverges at short 
distances the effects on the short distance wave function 
are not out of control; the relevant scale is provided by 
the binding energy of the bound state. 

Of course, while the outcoming decay ratios seem to 
reproduce the experimental numbers above the critical 
radius, it is not completely obvious whether fixing the 
short distance cut-off radius is fully justified, as the decay 
formula corresponds to the strict infinite mass and point- 
like limit of the related electromagnetic matrix element 
7 — > qq. Actually, the momentum dependence and re- 
coil effects due to the incoming virtual photon have been 
dealt with in Rcf. 47] where the corresponding modifi- 
cation to the standard decay formula decreases its value 
substantially. In appendix [B] we collect some relevant 
formulas. In our case we see that, as expected, there 
is an increasing suppression for excited states, although 
in the ratio this effect is negligible as can be seen from 
Table EH1 

In summary, the analysis shows that even if one can 
apply the renormalization to any long distance approxi- 
mation to the BSM potential, due to the existence of the 
scale ro, one should not take r c — > 0, but rather to keep 
Tc ~ fo to comply with the dynamics. 



B. Singular repulsive interactions 

To analyze a situation where renormalization condi- 
tions cannot be applied let us consider the case where in 
the power series expansion Eq. (|44p one term happened 
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State 


^-Wcisskopf 


^•Lcinweber 


J/i> 


1.00 


1.00 


i>{2S) 


0.58 


0.58 


i>(3S) 


0.41 


0.42 




0.31 


0.32 



TABLE VII: Relative n 3 Si — > e + e~ decay ratios with respect 
the ground state J/ip — ► e + e~ decay calculated through the 
van Royen-Weisskopf formula and taking into account mo- 
mentum dependence and recoil effects according to Ref. [47t ] . 



to be singular and repulsive at short distances, say for 
n = 3 (unlike the BSM discussed above) . Then the solu- 
tion would be of the form 



u(r)^r^ \c 1R e^r +C 2R e- 2 W 



(47) 



where we have a converging and a diverging exponential 
at the origin. 

Clearly, if we impose normalizability of the bound 
state, the regular solution must be chosen. In such a case 
the bound state energy is predicted as it usually happens 
for the standard non-singular quantum mechanical po- 
tentials. Conversely, if the bound state energy is chosen 
arbitrarily then the wave function is non-normalizable at 
short distances. In this sense, the rigorous results given 
by Eq. (|43j) . valid in the infinite mass limit, which im- 
ply that all inverse powers are attractive singularities are 
most welcome, because this would imply that the prob- 
lem of a repulsive singularity of a truncated potential is 
related to the truncation itself. As we have mentioned, 
these violations of the conditions featured by E q. (H3l are 
also observed in perturbation theory (see e.g. (15j|). 

Finally, it is not obvious what might be the situation 
for relativistic corrections which on purely dimensional 
grounds are of the general form 



a renormalization boundary condition scheme. This ap- 
proach avoids explicitly the introduction of phenomeno- 
logical form factors used in the past which are usually 
fine-tuned to low-lying mesonic observables. 

The renormalization viewpoint actually uses these low- 
lying properties as input parameters and sidesteps the 
fine tuning, as illustrated in several examples. For the 
phenomenologically successful model of Ref. [25| where 
ad hoc form factors are introduced as regulators, we find 
an almost perfect agreement with the renormalization ap- 
proach. This result provides confidence on the way the 
original model took into account the unknown short dis- 
tance dynamics. The use of the renormalization scheme 
allows us to further study the dependence of the states on 
the model parameters in a cleaner way since the regulator 
dependence has been removed. 

We find that the mass of the excited states strongly 
depends on the string tension parameter as one would 
naturally expect. Further, there is also a remarkable in- 
sensitivity to the strong coupling constant entering the 
one gluon exchange contribution of the potential. This 
avoids a great deal of unphysical fine tuning which sug- 
gested taking for this parameter unnaturally large values 
a s = 0.3 — 0.4. In contrast, the leptonic decay widths 
present a marked dependence on the strong coupling con- 
stant, as expected likewise because these are short range 
observables. Of course, a more complete description of 
charmonium systems should also take explicitly into ac- 
count the contribution to the open charm channels which 
will provide mass shifts and widths for unstable quark- 
antiquark bound states. In summary, the renormaliza- 
tion approach traces quite naturally and explicitly the 
long and short distance dependence of hadronic proper- 
ties. 



Appendix A: Renormalization conditions and Error 
estimates 



V(r) 



1 



P 



1 n 



(48) 



Of course, Weyl ordering should be applied to make sure 
that the operator is hermitian [Hj]. However, as shown 
in (4^ . IHoj one can eliminate the p dependence in terms 
of L dependence at the expense of reshuffling the 1/m 
expansion. Ultimately, we end up with singular prob- 
lems which generally depend on the angular momentum, 
showing that one must possibly use a different renormal- 
ization condition for each partial wave. However, note 
that for this to happen one must have a singular attrac- 
tive singularity. This interesting topic deserves further 
investigation. 



VII. CONCLUSIONS 

We have reanalyzed the calculation of the 1 — - char- 
monium spectrum in a constituent quark model using 



The boundary condition allows to connect the bound 
state to the excited states by the matching condition. 
Here we estimate what is the error in the excited mass 
as a function of the the cut-off radius. This means that 
we start with the ground state energy, Eq, which is fixed 
throughout and compute the logarithmic derivative at 
r c . Then, we perform a calculation of the bound state 
energy, E n {r c ), which obviously depends on the cut-off 
radius. Let us denote u n (r, r c ) the excited wave function 
where the dependence on r c is explicitly shown. Under 
an infinitesimal change, r c — > r c + Ar c we get 

- A< + UAu n = AE n u n + E n Au n , (Al) 

Note that the variation is defined for a fixed value of r, 

Au n (r, r c ) = u n (r, r c + Ar c ) - u n (r, r c ), (A2) 

and hence 

Au' n (r, r c ) = < (r, r c + Ar c ) - <(r, r c ) , (A3) 
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where here, the prime denotes derivative with respect to 
the r variable. Therefore, if we use the boundary condi- 
tion 



u' n {r c ,r c ) = L n (r c )u n (r Cl r c ), 



we get 



Au' n (r, r c ) 



- A(L n (r c )u n {r c ,r c )) 
9u' n (r, r c ) 



dr r 



(A4) 



(A5) 



Ar, 



From here we obtain 



\ L n + L n - U + E n )u n = -— h L„ — — , 



(A6) 



and multiplying this equation by u n and subtracting the 
original equation multiplied by Au n we get 



A<m„ + <Ait„ = A£ n ^, 



(A7) 



which integrating from r c and infinity and using that for 
a bound state u n (r, r c ) — > at large distances we get 



Ait' u n + u'Au„ = AE n 



(A8) 



On the other hand, since the ground state energy is 
fixed, AE = 0, we get 



— Au' uq + u' Auq = 0, 
and taking into account 

u n (r c ) 2 



AE n = (Li +LI-U + E. 



J r °° u n (r) 2 dr' 



(A9) 



(A10) 



so that using L n = Lq we arrive at 

A(E n -E a ) = (E n -E )-^^—Ar c . (All) 

J,, c u n (r) z dr 

Actually, integrating we get 

poo 

[E n (r c )-E ] dru n {rf = const. (A12) 



(A13) 



For a normalized state we have for small r c 
E n (r c ) -E = [E n - E ) [1 + r c u n (r c ) 2 



Therefore for a regular potential with a non-trivial 
boundary condition u n (0) ^ the error is at least lin- 
ear. For a singular and attractive potential, l/r n , the 
error is 0(rl +n ^ 2 ) up to some oscillations. As we see, 
the convergence is from above and proportional to the 
energy difference as well. This means that the effect of 
putting a finite cut-off fixing the ground state energy is 
repulsive and increases with the excitation energy. 



Appendix B: Momentum dependence and recoil 
effects as corrections to leptonic widths 

Momentum dependent effects could reveal significant 
corrections to the theoretical leptonic width. The ex- 
pression for e + e~ decay width of S-wave states in the 
center-of-mass frame of the meson and taking into ac- 
count those effects can be written as [13] 

r(e+e-) = 16a 2 e^^(3 £ ; e 2 -|fc| 2 )x 



E q + m q \p[ 

TO Eg 



3(E q + m q f 



W\)d\p\ 



(Bl) 

where e q is the quark charge in units of the charge of the 
electron, M is the mass of the meson and 



E„ 



E, 



^) 1/2 > 

1/2 



(B2) 



are the quark and lepton energies respectively. 

The static limit of Eq. (|Bip may be obtained by con- 
sidering the nonrelativistic limit 

\P 12 



< 1 

and in this case the integral in Eq. (|B1[) reduces to 
2 



(B3) 



m q M 2 Jo 



ip(\p\)W\ 2 d\p\ 



Recalling the Fourier transform 
1 



tp-r 



(2tt) 3 / 2 

and evaluating at the origin, we have 



ip{p)dp 



*=0) = 



(2 



(B4) 
(B5) 

l HW\)d\p\, 
(B6) 



so the integral above, Eq. (IB4I) . becomes 
2 



m a M 2 



i>(\p\)W\ 2 d\p\ = ^4^W= 0)| 2 . (B7) 



m q M 2 



r(e + e 



-|0(r = O)| 2 , (B8) 



With the kinematic relationships, the static limit of 
Eq. dHJ is 

_ 167ra 2 Q 2 

static M 2 

which is the well-known van Royen-Weisskopf formula. 

Figure [TT] shows the comparison of the relative n 3 Si — > 
e + e~ decay ratios with respect the leptonic width of the 
ground state, J /if) — > e+e~ . The ratios are calculated 
by van Royen-Weisskopf and Leinweber formulas. Left, 
middle and right panels of the figure show the results for 
the different terms contributing to the long distance ex- 
pansion of the bosonic string model truncated to a given 
order (LO, NLO and NNLO). 



17 




<c (fm) r c (fm) r c (fm) 



FIG. 11: Relative n 3 Si — > e + e decay ratios with respect the leptonic width of the ground state, J/ip — > e + e which are 
calculated by van Royen-Weisskopf (solid line) and Leinweber (dashed line) formulas. Left panel: LO. Middle panel: NLO. 
Right panel: NNLO. 



Appendix C: Short distance analysis of wave 
functions 

1. Single channel case 

In this appendix we summarize a few interesting facts 
concerning the short distance behaviour of wave func- 
tions. In the OGE case an attractive Coulomb like be- 
haviour holds, so that at short distances the reduced po- 
tential reads 



(CI) 



so that for the nS state we have short distance behaviour 



u n (r) -> A r , 



2R R & \RJ 



+ B n r, (C2) 



which is a linear combination of the regular wave func- 
tion and the singular one. Note that since the ground 
state has a given energy the irregular component does 
not vanish. On the other hand, the orthogonality condi- 
tion, Eq. (J7J, implies A n /B n = Aq/Bq so that 



Un(^c) An 
— > 

u (r c ) A ' 



(C3) 



which shows that the ratio between wave functions be- 
comes finite as the cut-off is removed, as can be seen at 
the right panel in Fig. [TO] 

For a power-like short distance singular potential we 
may keep the strongest singularity 



(CM) 



The solution of the Schrodinger equation requires the 
reduced potential 



U(r)=2vVir)^-±(£ Y 



(C5) 



where for convenience the variable R — (2fiC n ) 1 ^ 4 with 
length scale dimensions has been introduced. At short 
distances the reduced de Broglie wavelength is given by 



which fulfills 



d\(r) n ( r \ f 



or 1 



« i, 



(C6) 



(C7) 



for r <C R. In such a case the WKB method can be 
applied [5l| yielding 



u(r) 



«WKB(r) 

A 

WW 11 



dr^J-U(r) + if 



, (C8) 



where A and tp are undetermined amplitude and phase 
which may be obtained by matching to the exact solution 
in the region r ~ R. In the case of the singular potential 
given by Eq. (|C4|) we have for the m-state 



u m {r) -t A T 



sin 



_/r 
2 \R 



(C9) 



However, the orthogonality condition, Eq. (J7J, imply 
fm — ipo- Thus we obtain 



u m (r c ) 
Uo(r c ) 



A m 

AT 



(CIO) 



which shows that the ratio between wave functions be- 
comes finite as the cut-off is removed. For n = 3 this is 
seen at the right panel in Fig. [TOl 



2. Coupled channel case 

Here we undertake the short distance analysis of the 
spin-orbit and tensor interactions. At short distances one 
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may neglect all terms and just keep the 1/r 3 singular 
contribution yielding 



± R ^\ u(r)\ =Q 
-w"(r)l r 3 lM -20 I \ w (r)j 



-u"(r) 



(Cll) 



This system can be diagonalized by going to the ro- 
tated basis 




cos a sin a \ / u(r) 
— sin a cos ay \w(r) 



(C12) 



where the new functions satisfy 



(C13) 



and the Ri and —R 2 are the corresponding eigenvalues 



Ri = -^ + 2V3>0, 
-R 2 = ~-2Vz<0, 



(C14) 



and the mixing angle is 9 = 1.1 , a rather small value. 
At short distances the solutions of Eq. (|C13|) could be 
analyzed via the WKB method as we have done in the 
previous section, but for this case we can undertake the 
analysis in terms of Bessel functions, which short distance 
is analytically known. Actually, the solutions of — y"{x)— 
y(x)/x 3 = are 



3 



y/xJl(2/y/x) = ^008(71-/4 + 2/^) + 



y/xYi(2/y/x) = -^cos(7r/4-2/v^) 

\/7T 



(C15) 



whereas the solutions of —y"(x) + y(x)/x 3 = are 



Vi/i(2/>/5) 



1 



2V^ 



(C16) 



last equation also shows that generally solutions will di- 
verge as e 2 V Rl / r at the origin, hence preventing the 
bound state normalization condition, unless C\r = 0. 
In such a case the normalizable solution may be written 

as 



u(r) — s-cosaCU ( — ) sin 




sin a C 2 r \ I < 



(C18) 



w(r) — >-sinaCyi — I sin 
R2 



cos a C 2 r 1^-1 ' 

The three independent constants appearing here for 
the regular solution Ca, C 2 r and ip correspond to fix the 
energy Mj/^, the asymptotic D/S ratio 77, and the nor- 
malization condition. Note that when integrating from 
large distances to short distances with a given bound 
state energy and an arbitrary r\ we would always have 
a contribution from the exponentially diverging solution 
since C\r ^ 0. Thus, the condition C\r — predicts 
i] from the bound state energy. The foregoing analysis 
shows that for the 1/r 3 singularity appearing here the 
solution is ambiguous and the bound state energy for a 
given state has to be treated as an input. Of course, 
the orthogonality requirement between different states 
implies that if one state, (u n ,w n ), has a short distance 
phase ip n and another state, (u m ,w m ), has a short dis- 
tance phase (p m , one has 

/>oo 

= 2/i(M„ - M m ) / dr (u n u m + w n w m ) 
Jo 



[U n U m - U n U m + W n W m - W n W n 



-TT-C A,nC A.m siO-ifn 
n 2 



9! 



(C19) 



whence cp n = ip m is obtained. This shows that all states 
are linked through the spin-dependent splitting provided 
the ground state energy is given. 



All this amounts to write the solutions in the suitable 
form 



r 

R~i 



C\ R e 



ui(r) 



v 2 (r) -> Ca ( ) sin 



+2 



C 2 Re~ 



(C17) 



The four constants appearing here, Cm, C^r, C 'a and 
(p reflect that the total order of the system is four. The 



Appendix D: Fitting parameters of the renormalized 

model 



It is interesting to see what is the renormalized model 
sensitivity to the parameters. Of course, at this level one 
should envisage the possibly non-negligible contribution 
of the missing mesonic thresholds and the induced mass 
shift. To this end we fit the model parameters from the 
experimental masses of 1 cc states [38| and estimate 
their theoretical uncertainties as well as an educated mass 
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shift guess. The way to perform the first task is of course 
debatable. This requires some compromise as to what 
are the errors attached to the masses within the present 
framework. Naively one would just take the quoted PDG 
errors on the experimental masses. However, as we have 
mentioned above, the present model incorporates thresh- 
olds but no explicit coupling to meson-meson channels. 
This affects more significantly the higher excited states 
via a subthreshold induced mass shift, which we take as 
a systematic error of the model. 

We decide to make a x 2 fit using as fitting parameters 
the quark mass, m c , a c and fi c related with the confine- 
ment strength and the strong coupling constant, a s . We 
define the 



{M exp {i) - M teo {i,p)) 2 



(Dl) 



where p represents our model parameters and the experi- 



mental data arc taken as M. 



AL 



4,(2S) 



3686.093±0.143MeV, 



= 4039.6 ± 42.25 MeV and M, 



4361 ± 



37 MeV, where as explain previously the errors are taken 
as the half-width of the state. By minimizing the \ 2 
function wc obtain the theoretical uncertainties from the 



corresponding covariancc matrix at the minimum. The 
outcoming values for the parameters are 

m c = 1862 ± 12.6 MeV (0.68%), 

a c = 524 ±43 MeV (8.2%), 

[i c = 88 ±7.2 MeV (8.2%), 

a s = 0.41 ±0.14 (34%), (D2) 

where we put in parenthesis the relative uncertainties in 
percentage. Here we clearly see that the highest uncer- 
tainty corresponds to the value of the strong coupling 
constant. This rather large insensitivity to the otherwise 
too large a s is a rewarding feature of the renormaliza- 
tion approach. Actually, quite natural values of a s are 
obtained. The spectrum at the x 2 — minimum is given by 



1>(2S) 
4>(3S) 



3687 ±80 MeV, 
4108 ±79 MeV, 
4348 ± 80 MeV, 
4586 ±66 MeV, 



(D3) 
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